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In this paper we study the thermalization of a spatially homogeneous system in a 
strongly coupled CFT. The non-equilibrium initial state is created by switching on 
a relevant perturbation in the CFT vacuum during At > t > —At. Via AdS/CFT, 
the thermalization process corresponds to the gravitational collapse of a tachyonic 
scalar field (m 2 = —3) in the Poincare patch of AdS§. In the limit At < 
the thermalization time is found to be quantitatively the same as that of a non- 
equilibrium state created by a marginal perturbation discussed in Ref. 5|. In the 
case At > A, we also obtain double-collapse solutions but with a non-equilibrium 
intermediate state at t = 0. In all the cases our results show that the system 
thermalizes in a typical time tx — ^p^. Besides, a conserved energy- moment current 
in the bulk is found, which helps understand the qualitative difference of the collapse 
process m 
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I. INTRODUCTION 

In the AdS/CFT proposed by Maldacena[lj], a complete understanding of the M = 4 
SYM (CFT) requires the information of the whole infinite tower of massive Kaluza-Klein 
modes on AdS 5 x S 5 as well as the low energy supergravity on AdSz,]2\. Massless, massive 
and tachyonic fields in AdS& respectively correspond to marginal, irrelevant and relevant 
perturbations in the boundary CFT. Thermalization of the strongly coupled system in CFT 
may be studied holographically by reconstructing the bulk geometry from the back-reaction 
of those matter fields j^]. 



The thermalization of a spatially homogeneous system initially in a non-equi 
state created by a marginal perturbation has been studied both perturbatively 



ibrium 

4] and 



numerically jsj]. The perturbation is excited by a scalar source coupled to a marginal com- 
posite operator. On the AdS side, one needs to study the black hole formation from the 
gravitational collapse of a massless scalar field in the bulk. It was shown that such a strongly 
coup.ed system thennalizes in a typical tune t T ~ ^ with T being the then„al-e q uihbn Um 
temperature [5|. is the typical time for the scalar field to fall from the boundary to the 
interior of AdS§ close to the location of the apparent horizon. For the other aspects of the 
;hermalization of such a spatially homogeneous system, the interested reader is referred to 



6|-|21| for CFT and (22|, |23| for classical field theory. 

If the boundary source is switched on only in a short duration, one may expect that 
the thermalization time tr should not (dramatically) depend on whether it is coupled to a 
marginal or relevant operator. However, on the gravity side the thermalization processes in 
these two cases respectively correspond to the gravitational collapse of massless or tachyonic 
scalar fields. It is interesting to see under what conditions the mass of the scalar field does 
not affect significantly the thermalization time. 

In this paper we study the thermalization process of a spatially homogeneous system. The 
system is initially in a non-equilibrium state created by turning on a scalar source coupled 
to a relevant operator Oa in CFT. The energy is injected into the system according to [3] 

^ro ( o 4) = -(OA>0; O ), a) 



where (Oa) is the expectation value of Oa and the overdots denote time derivatives. After 
the source is switched off, the system will finally thermalize through self-interaction and the 
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information about the initial state will be totally lost. In AdSd+i, it corresponds to the 
gravitational collapse of a tachyonic scalar field of mass m. Stability requires not m 2 > 
but m 2 > — The conformal dimension A of Oa is related to the mass of the scalar field 

byf 



A(A -d)=m 2 , i.e., A = ^ + y^- + m 2 = ^ + n (2) 



with 



d 2 

n=\J— + m 2 . (3) 

The scalar source induces a tachyonic scalar wave falling from the boundary into the interior 
of AdSd+i- It satisfies the boundary condition 

0(t, W )~^- A (o) (t). (4) 

The induced wave will eventually collapse to form a black hole, which is the gravity dual 
of a static plasma. At the end, it will be completely hidden behind the apparent horizon. 
The thermalization process in CFT will be studied by investigating such a collapse process 
in AdSd+i- In this paper, d — 4 and the scalar source 0( O ) is chosen to take the form 

(O ) (t) = -e- at \ (5) 
a 

where e and two parameters and At characterizes the duration of the source 

being turned on. 

This paper is organized as follows. In Sec. Hi] we give the equations of motion for a scalar 
field of arbitrary mass coupled to gravity in the Poincare patch of AdS$. The numerical 
scheme in addition to initial conditions and boundary conditions is also discussed in this 
section. Then, we study the propagation of a scalar field of arbitrary mass in AdS$ with the 
back-reaction to the bulk geometry being ignored in Sec. IHIl In Sec. IIVt we investigate the 
details about the gravitational collapse of the tachyonic scalar field with m 2 = —3. In Sec. 
fVT. we briefly conclude. The near-boundary behavior of the massless scalar field is presented 
in Appendix |A] In Appendix [Bj a conserved energy-momentum current in global AdS^ is 
defined, which is useful for understanding some details about the collapse process of scalar 

nn 

fields discussed in [9|, 110] . 
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II. EINSTEIN-KLEIN-GORDON EQUATIONS 

A. Equations of motion 

On the AdS side, we need to solve the Einstein-Klein-Gordon equations 

1 



-g 



d a {V=gg ab d b (p) - m 2 <p = 0, (6) 



1 d(d-l) 

ttab — -^dab-ti ^2 9ab ~ J-ab, { < ) 

where T ab takes the form 

T ab = 28^8^ - g ab [{d4>f + m 2 2 ] . (8) 

In the following, we take L = 1 and d = 4. For the spatially homogeneous system on the 
boundary M 4 , we use the Schwarzschild coordinates of the form 

ds 2 = -i (-fe~ 2S dt 2 + f- l du 2 + dx 2 ) , (9) 

where / and b are functions of t and u only. In this coordinate system, one obtains from 
eqs. 06]) and ([7]) the following equations of motion 



P=(fe' 5 V)', (10b) 
/ = ^uf 2 e- 5 VP, (10c) 

S' = In (V 2 + P 2 ) , (lOd) 



f 2 
t = —u 

J 3 



f{V 2 + P 2 ) + ^ 



2 

2 



+ -(/-!), (10e) 
u 



m 
u 

where the derivatives with respect to t and u are denoted respectively by overdots and 
primes, P = <p f and V = f~ 1 e S (p. 

B. Initial and boundary conditions 

Before the scalar source is turned on, the scalar field is assumed to vanish in the bulk. 
In this case, at the initial time t = ti the metric functions and the scalar field are as follows 

/ = 1, S = 0, V = and P = 0. (11) 



The boundary condition for solving (llOdj) . (HOaj) and (llObj) are respectively given by 



5 

■0 



6(t,0)=0, (12) 
V(t, u min ) = -2tee~ at2 u 2 ~^ and P{t, oo) = as u min 0. (13) 

C. Scaling symmetry 

The scaling symmetry of the equations of motion in eq. (jlOj) allows us to obtain solutions 
from a known solution as follows [5] 

4>\(t, u) = 4>(Xt, Am), f\(t,u) = f(\t,\u), and 6\(t,u) = 8(Xt, Am), (14) 

where (0, /, 5) denote the solution with the boundary condition given by 

<p{t,u min ) = ^e~ at2 (u min ) 2 ~ n (15) 

and ((fix, fx, 8\) with A > denote the one with 

\4-n 

<fi(t, u mm ) = -e-< x ^ (Xu mm ) 2 - n = ^e-^u 2 ^. (16) 

Or, equivalently, 

(a, e) (A 2 a, A 4 ~"e) . (17) 
As a result, we only need to study the dependence of solutions either on a or e. 

D. Numerical scheme 



The equations of motion in eq. (UOl) can be solved by the numerical scheme described in 
detail in Ref. jlj]: Given /, 5, V and P at t — tj, we first calculate V, P and / at the next 



time step tj + i by solving (llOaj) . (llObj) and (llOcR . Then, 5 at £,- + i is obtained by solving eq. 



(llOdl) . Given the initial conditions in ( ITT]) , the bulk metric at late times can be calculated 
by repeating the above two steps. In our numerical calculation, we use a grid in the u 
direction with uniform spacing du. We find that a stable algorithm requires the time step 



1 In our numerical calculations, we have also tried providing the boundary conditions by including higher- 
order terms in the power series solutions of eq. (fT0|) in eqs. (|41[) . (|42|) and (|43|) but we only find negligible 
difference. 
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dt < O.ldu. Also, du should satisfy du < Q.lu min to ensure that the boundary conditions 
can be correctly provided. 

For arbitrary mass, we can not in general provide the boundary condition at u = due 
to our choice of the coordinates in eq. (see eq. (Tj0)). In this paper, the boundary 
condition for numerical simulations is provided at u — u m i n <C 1 in eq. (TT3l . We find that 
our numerical results are insensitive to the choice of u m i n as long as u m i n < 0.01. For all the 
numerical results in the following sections, we choose u min = 0.001. 



III. SCALAR FIELDS IN AdS 5 



In this section we ignore the back-reaction of the scalars to the bulk geometry. In this 
case, the AdS$ metric, with f = I and 5 = in eq. (jSJ), has a killing vector 



which allows us to define a conserved total energy M(t, u) as follows 

8 8 r°° 



'0 c u \ u 

loo ' 



where 



£ o = JL f v 2 + P 2 + %d> 2 ) , and, £ u 



4 j., 



(18) 



(20) 



(21) 



3tt 3 \ u 2 J ' 3u 3 3u 3 

In Sec. IIVt we shall show that in the case - 1 M(t, 0) is a very good approximation to 
the black hole mass resulted from the gravitational collapse of the scalar fields. 

In such a fixed gravitational background, the equation of motion of the scalar field in eq. 
dSJ) or eq. (jlOaO reduces to 



m 



+ —<P + - 

u z u 



3,/ 



0. 



The general solution to the above equation is given by 



Re u 2 j pe-^ [Cj(u)J n (uu) + C Y (u)Y n (uu)] , 

Z7T 



(22) 



(23) 



2 That is, we have 



0= /rfx- 5 ^Va(T ab 6)= [<Px r^t a Z b T°»* +f -,u a ^ 

J J J u u —oo J —oo 11 



(19) 



with u a = and t a = — 
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where J n and Y n are respectively the Bessel functions of the first and second kind. For the 
boundary condition in eq. (pO 



with 



C Y (co) 



Re u 2 I ^e- iuJt C Y (u)Y n (oou) 

Z7T 



2 ^w77j^(0)( w ) and 0(o)M = / dte lwt <P m (t) = —e ^. 



3 

a 2 



(24) 



(25) 



Taking the explicit form of 0m) in eq. (J5J), we have 

CV(w) = 



eo; n e "a /7r\f 



2 n r(n) Va 



(26) 



At time £ ^> At = regardless of its mass the group velocity of the scalar wave in eq. 
is always equal to the speed of light. This can be proven as follows: (f) is not significantly 
suppressed only in the region \t — u\ < At at late times. Using the asymptotic expansion of 
Y n (C) at |C| > A 



(-i)' r c ; 



Yn(C) 



K) 2 



r [sinC - (-l) n cosC] 



if C>o 



(-i) 



' [ (2 - (-1)" i) sin C - (i + (-1)™ 2) cos C] , if C < 



(27) 



with n = 21 + (— l) n mod 2 in eq. f )24|) . after some algebra we obtain 

t — u 



At 



'o 



U){t — U) u){t-u) 
—1) cos — ^ — - + sin 



At 



At 



(28) 



where 



, 3 n\ „ / 3 n 3 ~\ , / 1 rA _ / 1 n 1 ^ 9 

' 4 + 2 ) 1 (j + 2' 2 ) + F 4 + 2 'M 4 + 2 ' 2 ' " C 



with 



C n = 



n 5 



(29) 



(30) 



2 n T(n) 

Here, iF x is the Kummer confluent hypergeometric function. Therefore, at late times all the 
scalars propagate at the speed of light in the bulk. 



3 Here, we only discuss the cases with n being integral and m 2 > —4. 
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Near the boundary u = 0, the power series expansion of the scalar field has the form 



d-A 



HO) + 0(2) U 2 + ■■■ + 4>(2n-2)U 2n - 2 + <\> { <2n)U 2n + ^n)^™ log V? + • ■ ■ ) , (31) 



where 0(2m) f° r 2m < 2n and ^(JJn) are determined by the boundary source 0( O ). The 
expectation value of the dual operator Oa is given byj"" 



(Oa) = 2n0 2ri 



(32) 



up to terms contributing contact terms in the 2-point function. In the next section, we shall 
focus on the tachyonic scalar with m 2 = —3. In this case, A = 3 and n = 1. The power 
series expansion of in eq. (124"1) near the boundary takes the form 

,2 



Re „ / 



(o) 



1 - ^ (2 7£ - 1 - log4 + 2 log uuj) u 2 + ■ 



u 



( 



\ 



4 

V V>(2) 



0) lOg M + + 




u 2 + 



(33) 



where je — 0.5772 is the Euler-Mascheroni constant and 4>2 dd or 0^ en are respectively an 
odd or even function of t defined by 



0(2? 



doj sm(cut)cj 0(o) (oo), 

1 f du 



g» = -(2 7j5 -l-lo g 4) (o) 



4 / 2tt 



0(0) (w)o; 2 logw 2 



(34) 
(35) 



Here, we have assumed that 0( O )(o;) is an even function of u. 



IV. TACHYONIC SCALAR FIELDS COUPLED TO GRAVITY 

In this section, we reconstruct the bulk geometry from the gravitational collapse of a 
tachyonic scalar field with m 2 = —3. In eq. fl7J), we implicitly assume that the VEV of the 
dual operator scales as N 2 . In this case, the back-reaction of the scalar field to the bulk 
geometry can not be ignored and it eventually collapses to form a black hole in the bulk. 
Such a collapse process is dual to the thermalization process in the boundary CFT. We 
study both narrow (^ < 1) and broad(^ > 1) waves. 
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A. Energy conservation 

Eqs. fllOcj) and (11 Pep give the conservation of the total energy in the bulk 

8 8 f°° 

-M(t,u) = -j u du£°=£«\l, (37) 



where the energy density and the energy flux along the u direction are defined by 
/ (V 2 + P 2 ) + 



m 2 2 



, and, = - JL/V'PV = -g^/0'0- (38) 



They reduce to those in eq. (121]) if the back-reaction of the scalar field to the bulk geometry 
is ignored. After the source is switched off, M(t, 0) is the black hole mass M^. This can be 
easily seen by rewriting eq. (llOep in the form 

7-iY 



u 4 



£°. (39) 



Integrating over u, one gets from the above equation 

poo f — 1 

M(t, 0) | t>At = / du£° = - lim J -— = M bh = (ttT) 4 (40) 

by Birkhoff's theorem. Here, T is the thermal equilibrium temperature of the boundary 
CFT or, equivalently, the Hawking temperature of the black hole formed in the bulk. 

In the boundary CFT, eq. (1371) also gives the conservation of the boundary energy in 
eq. ([1]). Now let us take m 2 = —3. Near the boundary the power series solutions to the 
equations of motion in eq. ( flQj) have the form 

8 = \<t>l)U 2 + ■ • • , (41) 

/ = 1 + a (2 )U 2 + (6(4) logw + a (4) )w 4 H , (42) 

= u (0(o) + 0(2) w 2 + ^(2)U 2 log u 2 H ) , (43) 

where the coefficients are given by 

2,2 



0(2) = 5^(0) , 6 (4) = 5 (0(0)0(0) + 0(0) + 20(0) J , ^(2) = ^ ( 3 0(O) + 4 0(O) J ■ ( 44 ) 



Equivalently, one has 



^» + ^" = 0. (36) 



£ Q = (£°,£") is referred to as the conserved energy-momentum current in the bulk. 
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FIG. 1. A typical tachyonic narrow wave (e = 0.5 and a = 60): At = 0.13 ~ ^ 18 . This figure 
shows u~ 3 / 2 as a function of u at different times. At t < 1.0 = 7. 7 At, the fixed-background 
solution in eq. (|24| agrees every well with our numerical solution. Then, the fixed-background 
solution, which can be very well approximated by the asymptotic solution in eq. (|28p . fails to 
describe the collapse process later on. At t > = 6.25, the scalar field is (almost) completely 
hidden behind the apparent horizon. 

Here the boundary condition in eq. (j4j) has been used and all the other coefficients of the 
power series can be expressed in terms of 0(o) and 0(2). Among them the nonsingular terms 
of eq. ( 157)) in the limit u — > give 

2 



*W = g 



*0(o)0(o) + 60(o)0(2) + 30(0) (60(2) + 0(0) J , (45) 

which is the time-evolution equation of the boundary energy density 

f l I <• \ ■ 3 1 • 1 

e(t) = - dt {0 3 j 0o = 0(0)0(2) - |G(4) + 7^(0) + ■ztfay ^ 



B. Narrow waves) At = 4= < m )' - < 1 

From Fig. [T]one can understand how a typical narrow wave propagates in the bulk. This 
figure shows w~^0 as a function of u at different times. At t < 7. 7 At, the fixed-background 
solution in eq. ( 1241 is a good approximation of our numerical solution. Then, due to the 
back-reaction of the scalar to the bulk geometry, the group velocity of the wave starts to slow 
down while the fixed-background solution still propagates at the speed of light, which can 
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FIG. 2. The location of the apparent horizon uq = -4? f° r m2 = — 3. Here, we fix e = 0.5 and take 
uq as a function of a. The red solid curve (Analytic) shows the analytical result in eq. (j53[) . which 
agrees very well with our numerical results for all a > 1, or, equivalent, At < 

be clearly seen from the asymptotic solution in eq. (128]) . As a result, one can use the fixed- 
background solution to calculate the thermal equilibrium temperature T and to describe the 
details about how the energy is injected into the CFT vacuum during At >t > —At. 



1. The thermal- equilibrium temperature T 

Let us first calculate Mf,h = i^T) 4 , the black hole mass formed in the bulk. Taking e = ^ 
as a small parameter, we write /, 5 and in the formj^] 

/=l + e/ 1 + e 2 / 2 + 0(e 3 ), (47) 

S = e6 1 + e 2 6 {2 )+0(e 3 ), (48) 

= e0 1 + e 2 2 + O(e 3 ). (49) 

The equations of motion in eq. (jTUl) can be solved order by order in e. At the first order in 



A = OA = o, (50) 



and e0i is given by the fixed-background solution in eq. fl24j) . Then, the second order 
equation 

h = y0'i0i (51) 
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gives the black hole mas 



,:„*7a(°°)_ 



= -0( 4 )(oo) = - lim — = - lim — - / dty'xfa, (52) 

u^O -U 4 u^O dll 6 J_ 00 

where a( 4 )(oo) = — ^ and uq = ^ is the location of the apparent horizon. Therefore, 
M(t,0)\ t>At of the fixed-background solution in eq. ( 120]) is approximately the black hole 
mass in the bulk. Inserting ( )33l) into ( I52jl and only keeping the terms even in t, we get 



/CO / A \ o /-co 47Tf 2 f 2 

eft (4(0)^ + -0 ( o)0^ 3 y dtiwfffi = ^^~- ( 53 ) 

By comparison with our numerical results in Fig. |2} we find that eq (1531) is a valid approxi- 
mation for At < ^r. Using <p m e Q- (IA1 [) as e</>i for the massless scalar, we have (4 ] 

16VT r .. , , odd 47T6 2 _ , 



Mw, = -— y ^0(0)0(4) = -3- oc e 2 , (54) 

which also agrees with our numerical results in Ref. j^]. 

The above scaling behavior of is a consequence of the scaling symmetry. Under the 
scaling transformation in eq. (JT] 



M bh A 4 M 6ft . (55) 
Therefore, at 0(e 2 ) one should has, for arbitrary n, 



2 



according to the scaling transformation in eq. flTTj). This scaling behavior should be inde- 
pendent of the shape of the scalar source. 



2. Energy injection 

Using the (approximate) analytic solution above, we discuss the procedure to obtain the 
boundary energy from the numerical bulk solutions in this subsection. Even though we are 
mainly interested in the thermalization time tr and e (strongly) depends of the shape of the 
boundary source, the lesson we learn from such a simplified case may be useful for studying 



the more complicated systems, say, in 24j. 134 1 . 



5 Here, we take t — 00 but because of the energy conservation in eq. ([37| it is the black hole mass after the 
scalar source is switched off, that is, at time t > At. 



13 




FIG. 3. The boundary energy density e(t). Fig. (a) shows the fit results of e with u s = 0.4 and 
n s = 24,26 and 100. We find that for n s > 26 they do not depend on n s , which agree with the 
result (Analytic) in eq. ([57|) . In contrast, for u s = 0.5 the fit results always depend on n s as shown 
in Fig. (b). Here (o,e) = (60,0.5). 

The energy is injected into the vacuum in the boundary CFT according to eq. that 
is, t 

e(t) = - I dt (6 3 \ 0o = -2 / ^0(0)0(2)- (57) 

J oo J oo 

Insert into the above equation </>( 2 ) denned in eq. ( 133]) . and we find the following relation for 
different (a, e) 

e 2 1 • 

e(t) = -e(v^t)--0 2 o) loga (58) 
a 4 w 

with 

= l-^i ■ (59) 
The term oc log a is due to the conformal anomaly 26]. In order to get the boundary energy 



density one only needs to calculate (p^) from known solutions of eq. (fTUl) and then integrate 
over t in eq. fl57|) . 

Numerically, we find that it is easier to obtain e(t) in the following way: By inserting eq. 
( |33l) into eq. ( loTT) . near the boundary one has 

e 2 /2 = a(2)M 2 + (&(4) logw + a (4 ))w 4 H , (60) 

6 Note that the energy density may take negative values at some intermediate time, which can be understood 
from the second term on the right-hand side of eq. (|58|) . One may choose a renormalization scheme 
different from that used in Ref. [3| to make e be positive all the time. The interested reader is referred to 



25 1 for a recent discussion of holographic renormalization in the cases with m 2 = —3 and —4. 
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where the coefficients are given by 



0(2) = ~0(O), &(4) = 3 (0(0)^(0) + 0(0) j , ( 61 ) 



2 ^ h 2 

g0(o) s 0(4) = 3 

«(4) = J (0(0) + 40(0)0(2) + 8 / ^0(0)0(2) ] • (62) 



Accordingly, 



3 1 12 



e(t) = 0(0)0(2) - ^0( 4 ) + ^0(0), (63) 

which is eq. ( |4"6"]l up to 0(e 2 ). The coefficients am and 0( 2 ) can be obtained by the least- 
squared fit of the power series expansion of / and near the boundary to our numerical 
results. Here, the power series expansion of / and takes the form 

lis. 

2 

f s = l + a {2) u 2 + ^ ( b (2k) log u + a (2k)) u 2k , (64) 

(65) 



u 



k=2 
-1 

,2k 



0(0) + (^( 2fe ) logu2 + ^( 2fc )) ^ 

fe=l 

where n s > 4 is an even integer. The coefficients 0(o) , ip(2) , a(2) and 6( 4 ) take values according 
to eqs. (jSJ), ( |33l) and ( 16T|) but others are to be obtained from the least-squared fit to our 
numerical solutions at u — (u m i n ,u s ). u s > u m i n and n s are chosen such that the fit results 
of 0(4) and 0( 2 ) do not change by increasing n s . Then, the fit result of e is obtained by eq. 
(|63|) or eq. (|46l). 

Let us take for example the narrow wave with (a, e) = (60,0.5) (see Fig. [1]). In this 
case, we find that the fit results of e do not depend on u s in the range u s = (0.3,0.45). 
For u s = 0.3, 0.4 and 0.45, it respectively requires n s > 18, 26 and 42. Fig. E^a) shows our 
fit results with u s = 0.4, which agree with the analytic solution in eq. ( 1571) for n s > 26. 
However, for u s > 0.5 the fit results always depend on n s (see Fig. E(b)), which indicates 
that it = 0.5 should be outside of the radius of convergence of the metric functions. We also 
find that the fit result slightly depends on n s for smaller u s . It is due to the error of order 
u min = 10~ 6 introduced by the boundary condition in eq. ( 1T31) . 



3. The thermalization time 



The details about the collapse process of the tachyonic scalar is shown in Fig. [TJ Like the 
massless scalar Q], it shows that the system in the boundary CFT thermalizes in a top-down 
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FIG. 4. The thermalization time = %■ Here gt is taken as a function of At T. For the narrow 
waves with At < , the relative difference of gt for the tachyonic and massless scalars is less than 
0.8%. 

manner js]. Since the naive perturbative calculation breaks down at tx > t ^> At, we instead 
resort to numerical simulations to study the thermalization time tx- Here is denned to 
be the first time when 

min/(t T ,u) = 0.01, (66) 

which has an interpretation in terms of a spacelike Wilson loop (W(l ~ ^)) in the boundary 
CFT[5|. 

First, let us fix e = 0.5 and numerically calculate frp and T as functions of a = a^°' 5 \ which 
are denoted respectively by tj?' 5 ^ and T^ 5 ). Then, let us use the scaling transformation in 
eq. ( fT4"l) to get the result for arbitrary (a, e): 

(t T ,T,a,e) = ||^-,Ar( a5 ),A 2 a( a5 ),0.5A 4 -"j , (67) 

that is, the scaling parameter A is given by 

A = (2e)^ , (68) 

and the thermalization time £t is 

tr = % with * ^T<^ = ( 69 ) 
1 iru 



In order to compare with the results of the massless scalar in Ref. 



5] we take gt as a function 



of At T. The result for g t at 0.3 > At T > 0.015 is shown in Fig. H For At < we find 
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FIG. 5. Two-stage collapse in the bulk. These figures show how £° and / of a double-collapse 
solution with (a, e) = (0.3, 0.5) evolve over time. The qualitative features of the time evolution of 
those two quantities can be understood by the conservation of £ in eq. (|37p and the back-reaction 
off to / in eq. (I55]). 



that the relative difference of g% from that of the massless scalar is less than 0.8% and the 
mass term does not play an important role for very narrow waves. In summary, we conclude 
that in the cases with At < ^ the system in the boundary CFT thermalizes in a typical 
time tj> ~ h. 



C. Broad waves(At > ±): ^ > 1 

At the end of this section, we briefly discuss the double-collapse solutions of the tachyonic 
scalar field. It also undergoes two-stage collapse in the bulk, which is qualitatively the same 
as the massless scalar field in Ref. [5|: The energy injected into the bulk from the boundary 
before and after t = is respectively responsible for the two peaks in the bulk energy 
density S°. As a result, there are two local minima / ~ 0.01 respectively at u = ul and Mo- 
Fig. [5] shows how S° and / of a double-collapse solution with (a, e) = (0.3,0.5) evolve over 
time. At each time, / reaches its minima around the locations of the peaks in £ °. And the 




FIG. 6. Single-collapse (Fig. (a)) and double-collapse (Fig.(b)) solutions. They are respectively 
the analogs of the single-collapse and double-collapse solutions of the massless scalar fields in Ref. 



unsmoothness near its minima can be qualitatively understood by eq. f l3"9~j) : Near the peaks, 

E^=(L=1^L. (70) 

At t > At, the total bulk energy M(t, 0) = M bh is conserved according to eq. (|37|) . which 
governs the time evolution of S° shown in Fig. [51 

However, as shown in Figs. [5] and EJ / at t — is far from that of the AdS black hole 
metric near the boundary^]. This is in sharp contrast with the double-collapse solutions of 
the massless scalar field. Such a difference can be easily understood according to eqs. (jHj) 
and (1A7I) . which show that / ~ 1 + |0^-u 2 for m 2 = —3 while / ~ 1 — |0 2 o ^m 2 for m 2 = 0. 
At t — the system does not look like thermalized for either local or non-local operators. 

To conclude this subsection, we discuss the thermalization time for broad waves. In this 
case, we define the thermalization time £t by subtracting out At from that defined in eq. 
( 166L That is, we study the thermalization of the isolated system after the source is switched 
off. For the solution with (a, e) = (0.3,0.5), we find t T = 1.85 = ^ with u = 1.09. We 
find that this is also parametrically true for other double-collapse solutions. Therefore, we 
conclude that the typical thermalization time for such a strongly coupled system is always 



7 For Umin = 0.0001/0.001, we find the fit result of e always depends on n s and u s . It indicates that even 
u s ~ 0.1 is outside of the radius of convergence of the metric functions. The calculation with a smaller 
u-min is numerically too demanding for us. 
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of order ^ no matter whether the non-equilibrium initial state is produced by marginal or 
relevant perturbations. 

V. CONCLUSIONS 

In this paper, we study the thermalization process of a strongly coupled system via 
AdS/CFT. The non-equilibrium initial state of the system is created by turing on a scalar 
source during At > t > —At, which is coupled to a relevant composite operator. On the 
gravity side, we study the gravitational collapse of a tachyonic scalar field with m 2 = —3 
in AdS$. In the case with At < ^ we find that the black hole mass My, = (7rT) 4 oc (-^) 2 
with e = - being the amplitude of the scalar source. In this case, the injection of energy 
into the CFT vacuum can be very well described by perturbative calculations. Moreover, for 
At < ^P, the thermalization time tr is found to be quantitatively the same as that of the 
non-equilibrium state created by the marginal perturbation in Ref. |5|. For the case with 
At > -j, a non-equilibrium intermediate state at t ~ is found in double-collapse solutions, 
which is very different from that of the massless scalar field. In all the cases we find that the 
system, after the scalar source is turned off, thermalizes in a typical time^ _~ rr^- Such a 



rapid thermalization seems typical of such a strongly coupled CFT (see 27H4l| for various 
other non-equilibrium initial states). 

After the source is switched off, one has 

£ u (t,0) = = S u {t, oo). (71) 

Both the boundary at u — and the hypersurface at u = oo act as reflecting mirrors to 
the scalar field. As a result, its total energy M(t, 0) is conserved. In our case, as shown in 
Fig. [TJ the scalar field always propagates inwards and eventually hides behind the apparent 
horizon. It can not reach u = oo. In contrast, in the global coordinate the scalar field can 
travel from the boundary/center to the center/boundary in finite time. In this case, both 
the center and the boundary act as reflecting mirrors. The scalar field may collapse to form 
a black hole after oscillating between the boundary and the center several times joj. The 
conserved energy- momentum current S a = (S°,S X ) in global AdS± is defined in eq. (1B3I) . 
which is useful for understanding the qualitative features of such a collapse process (see Fig. 
[7]). A more comprehensive description in terms of S a would help understand better such 
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collapse processes in global AdS. Besides, many discussions in Sees. |TT] and [TTT] are also 
applicable to massive scalar fields (with m 2 > 0), which are relevant for the thermalization 
of the non-equilibrium state created by irrelevant perturbations. However, is divergent 
near the boundary u = (see eq. flU)) and the back-reaction of the matter field to the bulk 
geometry requires a more careful analysis[2|, y(. Moreover, in our double-collapse solutions, 
at t = the speed of light-like geodesies ^ = fe~ s ~ 1 + \<j) 2 ^u 2 > 1 at u ~ 0. It 
would be interesting to learn more details about such an intermediate state as well as the 
thermalization process by studying other probes such as heavy quarks [421]. dileptons[15| or 
prompt photons [16j. We leave those unanswered questions for future research. 
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Appendix A: Near-boundary behavior of the massless scalar field 



For the massless scalar field, A = 4 and n — 2. The power series expansion of the 
fixed-background solution in eq. ( 1241) near the boundary is as follows 

du 



where 



Re / — V(o)(w 



1 + 



U 2 U 2 U! 4 U 4 



U 



0(0) - -J0(O) + 



4 64 
/ \ 

<t>l% n + 0ffi-H> W log u 2 

\ *(4) / 



( 3 -4 7s -41og^) + 



(0) .odd 



"32"' ^ 

U 4 ) 



du 
32tt 



sin(o;t)a; 0( O )( W ) 



ieven (0) 



(3 _4 7B + 41og 2)-l J ^e-^V (0) (a;)u; 4 loga; 2 . 



(Al) 



(A2) 
(A3) 



Near the boundary, the power series solutions to the equations of motion in eq. f lTUj) take 
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the form 8 



3 /(0) ' 72 
where the coefficients are given by 

1 ■ 



f = 1 + a(2)U 2 + (5(4) \ogu + 0( 4 ))w 4 + • • • , 

+ 0(2) U 2 + 0(4) M 4 + "0(4) M 4 logW 2 + ■ ■ 
1 

' "'^(0) "I" d V(0) u ^(U)r(o) / • 



6 = tU 2 0(o) + ^ 4 ( 16 0fo) + 3 0(o) - 6 0(o)0 (3; 



1(2) 



3^0)1 6 (4) 
1 •; 



6 V 4 0?O) + V,0:, 



0)0(0) 



(2) = -^0(0), ^(4) = 



l(4) 



32 V r (°) 

The nonsingular terms of eq. (l3~Tj) in the limit k -)• gives 



04 



As a result, we have 



1 

36 



1920(4)0(0) + 320^ O )0(o) + 60(o)0(oj ~~ 3 ^(o)0(J 



e(t) = - / dt(6. 



-4 



^0(0)0(4) 



7«(4) + 7^0(0) - ^0(0) 



t(3) 



16 



'(0) 1 6 



'(o)- 



(A4) 
(A5) 
(A6) 

(A7) 
(A8) 

(A9) 



(A10) 



Appendix B: Gravitational collapse of massless scalars in global AdS± 

The metric of AdS± in the global coordinate takes the form 9|] 

1 



ds' 



COS 2 X 



(-fe' 25 dt 2 + f^dx 2 + sin 2 xdn 2 ) , 



(Bl) 



where the metric functions / and 5 are functions of t and x. Here, the boundary of AdS± 
locates at x = ~. It takes only t = ~ for light to travel from the boundary/center to the 
center/boundary. From eqs. ([6]) and (j7j), one has 



^M(t, n, x 2 ) = ^ f 2 £° = £*{t, Xl ) - £ x (t, x 2 ) 



V = 



1 



tan 2 x 



(tan 2 xfe- s P)' , P = (fe^V)' , 



5' = - sin x cos x (V 2 + P 2 ) J = ^^-£*, 

sm x 



(B2a) 
(B2b) 
(B2c) 



At t = 0, / = 1 + ( log U + 0( 4 ) J M 



For the double-collapse solutions, one always has \a&\ 3> 



2f logu 



at ul > u > u m in- This makes / look like that of the AdS black hole at t = 0. 
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0.6 




(a) (b) 



FIG. 7. Gravitational collapse of a massless scalar in global AdS^. Here, we show the time evolution 
of £ (Fig. (a)) and £ x (Fig. (b)) of a solution which collapses to form a black hole at t ~ it. At 
t = 0, the matter field locates mostly near the center (x = 0). Then, the matter field propagates 
upwards as a result from the reflection by the center. Afterwards, it keeps going up before reaching 
the boundary, which forces it to fall inwards. At t ~ it, £° is more sharply peaked near x = 
than that at t = 0. The force from the reflection can not resist the gravity and the scalar field will 
eventually collapse to form a black hole. 

where the derivatives with respect to t and x are denoted respectively by overdots and 
primes, P = </>', V = / _1 e <5 </) and the conserved energy-momentum current £ Q = (^°,£ a; ) is 
denned by 

S° = tan 2 xf (V 2 + P 2 ) , and, £ x = -2 tan 2 xf 2 e~ s PV. (B3) 

After the boundary source is switched off, one has 

P(t, 0) = 0, and V (t, |) = 0. (B4) 

As a result, the total energy M(t, 0, |) is a constant according to eq. (lB2ap . £° and £ x help 
understand some details on how the scalar field evolves over time in the bulk. Fig. [7] shows 
the time evolution of S° and £ x of a solution with the initial condition given by P(0, x) = 
and V(0,x) = ^pexp ^— 1024 ^ n2x y it illustrates the general features about such a collapse 
process. All the discussions here can be easily generalized to the cases of massive/tachyonic 
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scalar fields in global AdSd+i or in Minkowski space (by taking the limit x — > 0). 
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